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Abstract 

Generalizing the well-known lilypond model ([U \TT\ [2]) we introduce a growth- 
maximal hard-core model based on a space-time point process of convex particles. 
Using a purely deterministic algorithm we prove under fairly general assumptions 
that the model exists and is uniquely determined by Under an additional station- 
arity assumption we show that the model does not percolate. Our model generalizes 
the lilypond model considerably even if all grains are born at the same time. In 
that case and under a Poisson assumption we prove a central limit theorem in a 
large volume scenario. 
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1 Introduction 

CO 

We consider a point process \I> = {(X n , T n , Z n ) : n > 1} on R x R + x /C , where 
R + := [0, oo) and K d denotes the space of all convex bodies containing the origin 6 R d 
in its interior. We interpret X n as the position of a grain (particle) with shape Z n that is 
born at time T n . Without interaction the n-th grain starts growing at time T n to form a 
grain X n + (t — T n )Z n := {X n + (t — T n )x : x G Z n } at time t > T n . It ceases its growth as 
soon as it encounters any other particle. This interaction can occur in two ways. First, X n 
can be covered by another grain by time T n , so that no growth can start at all. Second, the 
growing grain may touch another grain (which may itself be either growing or have ceased 
growing at an earlier time). As a result of this growth process a growth time R n e M + 
is attached to the n-th particle. The full-grown grain is given by Z* := X n + R n Z n if 
R n > 0. We use this notation also in case R n = 0. (Then Z* = {X n } has empty interior.) 
The growth protocol (informally) described above entails the following three properties. 

(i) The interiors of Z^ and Z* do not intersect whenever m ^ n. 
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(ii) If R m = then there is some n 7^ m such that R n > 0, X m e and T m > 
T n + inf {r > : X m e X n + rY n }. 

(iii) If R m > then there is some n 7^ m such that i? n > 0, T m + i? m > T n + R n and 
^ n z* ± 0. 

Property (i) says that {Z* : i? n > 0} is a hard-core system, while (ii) means that the 
growth of the n-th particle is inhibited by some other grain. Property (iii) says that 
any grain with positive growth time is stopped by some other grain that reached its 
final size at the same time or earlier. We summarize properties (ii) and (iii) by calling 
(X n ,T n , Z n ) (or Z*) an earlier neighbour of (X m ,T m , Z m ) (resp. Z^). Any grain must 
have at least one earlier neighbour. If {R n : n > 1} is a family of M + -valued random 
variables depending measurably on \& and such that (i)-(iii) hold almost surely, then we 
call ty* = {(X n ,T n , Z n , R n ) : n > 1} a growth-maximal hard-core model (based on \I/). 
Under certain assumptions on \1/ we will prove in Section [2] that the model exists and is 
uniquely determined by 




Figure 1: Two models based on the same configuration of positions but with different 
birth time scenarios. 

The lilypond model [U [HI [2] is a well-known special case of a growth-maximal hard- 
core model. In this case (T„, Z n ) = (0, B d ) for all n > 1, where B d is the (closed) unit 
ball in M. d . In the physics literature this was called touch- and- stop model, see pQ. In [5] 
the model was generalized to the case (T n , Z n ) = (0, B) for a general symmetric star-body 
B C M. d . Neither the case of random (possibly different) shapes Z n nor the case of a 
space-time driving process $ has been studied before. 

Figure [1] gives an impression of how the introduction of starting times alters the 
structure of the model. The set of grains is the same in both pictures. But on the 
left-hand side all birth times are zero while on the right-hand side the birth times are 
independent and uniformly distributed on the interval [0, 10]. We see that in the presence 
of birth times some germs might not be born at all. This in turn implies that the germs 
which are born have more space to grow. Figure [2] compares the length and cluster size 
distributions of the preceding two birth time scenarios. The distributions on the right- 
hand side have considerably higher probabilities at large values. 

The paper is organized as follows. In Section |2] we use an exteneded version of an 
algorithm in [2] to prove under fairly general assumptions that the model exists and is 
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Figure 2: Comparison of length and cluster distributions in different birth time scenarios. 

uniquely determined by \17 in a translation invariant way. In Section [3] we generalize the 
results in |U |2] and prove under appropriate stationarity and moment assumptions that 
there is no percolation in the model. In Section H] we consider an independently marked 
stationary Poisson process with constant birth times. Using stabilization arguments as in 
[9], we prove that the growth times R n , n 6N, satisfy a central limit theorem. 

2 Existence and uniqueness 

Consider \17 = {(X n ,T n , Z n ) : n > 1}, as introduced above, defined on an abstract 
probability space (Q, J 7 , P). Assume that $ := {X n : n > 1} is locally finite, that is 
:= card{m > 1 : X m G B} is P-a.s. finite for compact sets B C M. d . For any k > 1, 
let denote the k-th factorial moment measure of \17, the measure on (M. d x R + x K, d ) k 
defined by 

^ k \C) :=EY,*H(x 1 ,t 1 ,K 1 ,...,x k ,t k ,K k ) GCn$ fc } (2.1) 

for measurable sets C C (M d x M + x K, d ) k , where the * indicates that the summation 
is over all fc-tuples ((a?i, ti, Ki), . . . , (x k , t k , K k )) with pairwise different entries. We shall 
assume that 

„(*) < a k (\ d ®Q) k , k>l, (2.2) 
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where a > 0, X d denotes the Lebesgue measure on IR d and Q a probability measure on 
R + x K, d . We shall further suppose that the measure Q satisfies the integrability condition 

J (r 2d + l)(p(K) 2d + l)Q(d(r, K)) < oo, (2.3) 

where p(K) is the radius of the smallest ball circumscribing K G )C d . 

In this section we will prove the following existence and uniqueness result. 

Theorem 2.1. Assume that (12. 2p and (12. 3p ZioZd Then there exists a unique growth- 
maximal hard-core model = {(X n , T n , Z n , R n ) : n > 1} based on 

Under the given assumptions the theorem says that there is a family {R n : n > 1} 
of ^-measurable random variables such that \&* satisfies conditions (i)-(iii) given in the 
introduction almost surely. Moreover, if {R' n : n > 1} is another such family, then 
R n = R' n P-a.s. for any n > 1. 

Before proving Theorem 12.11 we provide a short discussion of assumption (I2.2p . For 
> 1 let denote the k-th. factorial moment measure of $ := {X n : n > 1}. If \I/ is an 
independent marking of $ with mark distribution Q, then 

z/ (fc) (rf(x 1? t la A 1; . . . , x k , t k , K k )) = /i (fc) (d(x 1; . . . , x k ))Q k {d{t x , K h ..., t k , K k )) 

and (12. 2 p is implied by /i^ fc ^ < a k (\ d ) k , k > 1. The most fundamental example of a point 
process $ with this property is a Poisson process with a bounded intensity function (e.g. a 
stationary Poisson process). Other examples are provided by certain Cox, Poisson cluster 
and Gibbs processes, see [2] for more details. 

Our proof of Theorem 12 .11 is based on a purely deterministic result inspired by Proposi- 
tion 2.8 in [2]. We start by introducing some notation. Define N as the set of all countable 
sets if) C M. d x M + x JC d such that card{t> G ip ■ v G B x W + x fC d } < oo for all bounded sets 
Bcl d and for all (x, s, K), (y, t,L) G if) we have (x, s, K) ^ (y, t, L) only if x ^ y. Any 
if) G N is identified with the (counting) measure card{t> G if) : v G ■} on M. d x M + x /C d . A 
point (x, s, i^) G ■?/> is interpreted as a groin with shape K located at x and born at time 
s. As usual we equip N with the smallest a- field making the mappings if) i— > if)(C) 
measurable for all measurable C C M. d x M + x JC d . The point process \I/ is then formally 
defined as a random element in N. 

Consider if) = {(x n ,t n , K n ) : n > 1} G N. We call a subset ■?/>* = {(x n ,t n , K n ,r n ) : 
n > 1} of M. d x M + x JC d x M + with projection ^ on M d x M + x /C d a growth-maximal 
hard-core model based on if) if the interiors of x rn + r m K m and x n + r^A^ do not intersect 
whenever m ^ n and any point in if)* has at least one earlier neighbour. If r m = there 
is some n ^ m such that r n > 0, x m G x n + r n A n and t m > t n + inf{r > : x m G 
x n + rK n }. If r m > there is some n ^ m such that r n > 0, t m + r m > t n + r n and 
(a; m + r m K m ) n (x n + r n K n ) ^ 0. For an earlier neighbour (x n , t n , K n ) of (x m , t m , K m ) we 
will equivalently use the term stopping neighbour because of the geometric interpretation of 
our model given in the introduction. Furthermore we will refer to the function R : if) — > R + 
defined by R(x n ,t n , K n ) := r n , n G N, as a hard-core function on if) if if)* is a hard-core 
model. If if)* is moreover growth maximal we call R a growth-maximal hard-core function. 
If R is a hard-core function on if) we call two different points u = (x, s, K),v = (y, t,L) G if) 
neighbours w.r.t. R if (x + R{u)K) D (y + R(v)L) ^ 0. 
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For convex and non-empty subsets K, K', L, L' of ~R d we define 

a(K, K', L, Lf) := inf {r > : (K + rK') n (L + rL') ^ 0} (2.4) 

and use abbreviations of the type a(x, K,y, L) := a({x},K,{y},L), x, y e M d . For two 
grains u := (x,s,K) and v := (y,t,L) we define the first contact time d(u,v) as follows. 
If s < t we set 

d(u, v) := l{a(x, K,y,0) <t — s}(s + a(x, K, y, 0)) 

+ l{a(x, K,y,0)>t- s}(t + a(x + (t - s)K, K, y, L)). (2.5) 

In the case covered by the first summand, the growing grain u reaches the point y not later 
than the birth time t of grain v. In the second case both grains start to grow and meet 
at time d(u, v). If s > t we define d(u, v) := d(v , u), so that d(-, •) becomes symmetric. In 
the following and later we write a A b := min{a, b} and a V b := max{a, 6} for a,i»6l. 

Lemma 2.2. Le£ R be a hard-core function on ip £ N. In addition suppose u = (x, s, K) 
and v = (y,t,L) to be neighbours in ip w.r.t. R and R(u) A R{v) > . Then 

(R{u) + s) A (R(v) +t)< d{u, v) < (R(u) + s) V (R(v) + *). 

Here the second inequality is strict, provided that R(u) + s ^ R{v) + 1. 

PROOF: We assume w.l.o.g. s < t and prove the assertion by contradiction. First 
suppose that (R(u) + s) A (R(v) + 1) > d(u, v). If a(x, K,y,0) <t — s then 

R(u) + s > d{u, v) = s + a(x, K, y, 0), 

where we have used (12. 5p to get the identity. Since y e x + a(x, K, y, 0)K, this contradicts 
R(v) > and the hard-core property. If a(x, K, y, 0) > t — s, then d(u, v) = t + r, where 
r := a(x + (t — s)K,K,y,L). Therefore R{v) > r. But the definition of a(-) implies 
(x + (t — s + r)K) n (y + rL) ^ and hence R(u) < t — s + r = d{u, v) — s, by the hard-core 
property. This contradicts our assumption R(u) + s > d(u,v). 

Now suppose (R(u) + s) V (R(v ) +t)< d(u, v). If a(x, K, y, 0) < t — s, then 

R{y ) +t< d(u, v ) = s + a(x, K, y, 0) < t, 

contradicting R{y ) > 0. If a(x, K, y, 0) > t — s we obtain 

R(v) + t < d(u, v)=t + a(x + (t- s)K, K,y,L) = t + r, 

so that the definition of a(-) implies (x + (t — s + r)K) PI (y + R{v)L) = 0. Since u and 
v are neighbours we get R(u) > t — s + r = d(u, v) — s, contradicting our assumption 
R(u) + s < d(u, v). 

To prove the final assertion we again argue by contradiction. As the case a(x, K, y, 0) < 
t — s can be treated as above we assume that a(x,K, y,0) > t — s. Assume first that 
d(u, v) > R{u) + s > R(v) + 1. As before we then get R(u) > d(u, v ) — s, a contradiction. 
Let us finally assume that d(u,v) > R(v) +t> R(u) + s. Then the definition of a(-, •) 
implies that R(u) > t — s + r = d{u, v) — s, contradicting d{u, v) > R(u) + s. □ 
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The translation of if) G N by y G M d is defined by y+if) '■— {(y+x, s, K) : (x, s, K) G ?/>}. 
A set H C N is translation invariant if y + if> E H for all if> G and all y G M d . For a 
mapping R : N x R d x M+ x /C d M+ and ^ G N we define Cl d xK + xlC d xl + 
by 

flt(^) := t ; K, R(iP, x, s, K)) : (x, s, K) G if)} 

If i? is translation invariant in the sense that 

R{y + if),x + y,t,K) = R(if), x,s,K), if) G N, x, y G M d , (s, AT) G M+ x /C d , (2.6) 

then W(y + if)) = # f 0) for all (if), y) G N x ]R d . 

Theorem 12.11 is implied by the following Proposition and Lemma 12.41 

Proposition 2.3. There exist a measurable mapping R : N x M. d x R + x /C d — > R + and a 

sei H G A/" snc/i t/iat if if) E H then R?(if)) is the unique growth-maximal hard-core model 
based on if). Moreover H and R are translation invariant. 

The next lemma shows that H is reasonably big. 

Lemma 2.4. Let H be the set from Proposition \2.c\ Then P(\l/ G H) = 1. 

Proof of Proposition I2.3t Our strategy for the proof is as follows. We first 
construct H and R via a recursive algorithm. Growth-maximality and the uniqueness 
property will be proved in two separate lemmas. 

Take if> G N. Starting with if> := we will recursively construct an increasing sequence 
ipi C if), i G N, together with a sequence of hard-core functions i?j : ^ — > R + such that 
-Rj + i extends i?, for any 2 > 0. We will further define R(ip,u) := Ri(u) whenever u G ipi 
for some i. Afterwards this definition will be extended to all other u G if) in a suitable 
way. 

Let z G N. For u E if)\if>i we define 

z/j(u) := {m} U {t) G t/i \ (-?/>j U {m}) : v and w are mutual nearest neighbours 

in if) \ ipi with respect to d} 

We assume from now on that if) satisfies 

d(u,v) d(u' ,v'), u, u', v , v ' G if), {u, v} ^ {u', v'}. (2.7) 

Then Ui(u) contains at most two points. Furthermore let 

(Pi := {u G ip\if)i : n{u) {u}} (2.8) 

be the set of all mutual nearest neighbours in if) \ ipi (w.r.t. d). In the case (pi = we 
define if>%+x = if>i i.e. ipj = ipi, j > i. Otherwise we obtain if>i+x based on ipi in the following 
way. Set f(u,v) := inf{r > : (x + rK) H (y + Ri(v)L) ^ 0} + s for u = (x, s,K) Eip\ipi 
and v = (y, t, L) G ipi. For u G (fi define di(u) := inf^. ,_R i ( 1) )>o f(u, v) as well as 

s[(u) := d{u,v), v G z/ n (u) \ {u} and s"(w) := min d^v). 

vevi(u) 
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Now we add points of <pi to if>i by the following procedure to obtain ipi + i. 

Consider u = (x, s, K) G tpi so that there exists v = (y, t, L) G Ui(u)\{u}. Assume first 
that s-(w) < s'l{u). If s <t and a(x,K,y, {0}) < t — s we set Ri + i(v) : = and add v to 

If s < t and a(x, K, y, {0}) > t— s we add w and v to ^ and define Ri + \{u) := s^(rt) — s 
and := s^(u) — t. In case i < s we apply the same rules as before with the roles 

of u and v interchanged. Now we assume s'-(u) < s-(u). For all w := (z,r,M) with 
w G Ui{u) and c?i(w) = s"(u), we set i2j + i(ty) := s"(u) — r and add those points to fy. 
Applying this algorithm to every u G (pi yields the (i + l)-th hard-core function Ri+\ on 

We define 

oo oo 

^oo:=(J^ and V : = ^ W'oo = (2-9) 

i=l i=l 

Set R(ip,u) := -Rj(w) whenever u E ipi f° r some z. If card(y9 = 1 and {f G : R(v) > 
0} ^ we define R(ip, u) := inf„ e ^ ooiRWil , )>0 /(«, u) - s for u = (x, s, K) G 9?. 

Let H denote the set of all ip satisfying (12. 7p and for which the above algorithm yields 
a set t/'oo with card(^ \ ^00) < 1 and {v G ^00 : R(4>, v) > 0} 7^ 0. For all ^ e N \ if, we 
set ^(V 1 , •) = 0. Then R is a hard-core function. By construction the set H is translation 
invariant, and R is translation- invariant on H. Thus R is translation-invariant. We need 
to show that R(ip, •) is for any if) G H the unique growth-maximal hard-core function on 
if). To aid reading, we do this in two separate Lemmas. 

Lemma 2.5. Assume that if) G H. Then R := R(if>, •) is a growth-maximal function. 

PROOF: Let u = (x,s,K) G if). To show that u has an earlier neighbour. We 
distinguish the cases R{u) = and R(u) > 0. 

1. Assume that R(u) = 0. By the algorithm (introduced in the above first part of 
the proof of Proposition 12. 3p there is an u G ^00, that is u G ipi + i \ ipi for some % > 0. 
Moreover, the algorithm implies that there exists a v — (y, t, L) G if) \ {u} such that 
{u, v} = Vi(u) and 

<{v) < s'!(v). (2.10) 

We show, that v is an earlier neighbour of u. First we assume that v G if>j+i \ tfj f° r some 
j > i (the other case being v G if) \ ipoo)- Thus there isaw = (z, r, M) & ip\ (ipj U {u, v}) 
with {v,w} = Vj{v). 

Assume now that s'j(v) < s"(v). Since u, v are mutual nearest neighbours in if) \ ipi D 
if)\if>i and (12. 7p holds (since if) G H), we have 

d{v,w) > d{v,u). (2.11) 

We now claim R(v) > 0. Assuming on the contrary that R(v) = we obtain from 
R{u) = and (J23]) that t < s, 

a(x, 0, y, L) =< s - 1 (2.12) 

and 

t + a(x, 0, y, L) = d(u, v). (2-13) 
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Using the same argument for the pair (v, w) we obtain 

t + a(x, 0, y, L) = d(u, v) < d(w, v ) = r + a(y, 0, z, M) 

and, moreover, a(y, 0, z, M) < t — r. It follows that t + a(x, 0, y, L) < t, a contradiction. 
Recalling s'Jv) < s"(v) we derive from R(v) > the first equality in 

R{y ) = d(v, w) — t > d(v, u) — t = a(x, 0, y, L) 

and therefore 

x G y + R(v)L. (2.14) 

Together with (I2.12p this implies that v is indeed an earlier neighbour of u. 

Assume now that s'j(v) > s'Uv). By the algorithm there exists a neighbour w' : = 
(z', r', M') G iJJj of v with R(w') > 0. Further there are k < j and w" = (z", r", M"), such 
that w' and w" are mutual nearest neighbours in ip \ipk- Hence the algorithm implies the 
first identity in 

R( w >) + r ' = s' k (w') A 4V) < 4K) = d {w\ w") < d(w', v), 

where the final inequality comes from the fact that w', w" are mutual nearest neighbours 
in ip \ tpk D ip \ ipj . Since w' and v are neighbours, Lemma [2.21 yields 

R(v) +t > d{v,w'). (2.15) 

Again we have to distinguish two subcases. If k > i we get from the fact that u and v 
are mutual nearest neighbours in ip \ ipi D ip \ ipu that d(w',v) > d(u,v). Together with 
fl2TT5|) and fr2TT3|) this yields 

R(v) > d(v, w') — t > d(u, v) — t = a(x, 0, y, L). 

Hence (I2.14p follows. Assume now that k < i. To show that 

R(v) = a (y, L, z' + R{w')M', 0) > a(x, 0, y, L), 

we assume on the contrary that the inequality fails. But then the algorithm implies 
s'^v) > s'-(v), contradicting f)2.10p . 

Now we assume that \ ipoo . To show that v is an earlier neighbour of u we can 

proceed analogously to the case s'^{v) > s"{y) from above. 

2. Assume that R{u) > 0. We have to consider two subcases. If s'^u) < s"{u) 
there exists v = (y,t,L) G i/) such that {u,v} = z/j(w), R{u) = d(u,v) — s > and 
R(v) = d(u,v) — t > 0. Therefore R(u) + s = R(v) + t and v is an earlier neighbour 
of u. If s'^u) > s"{u) then % > 1 and there exists a neighbour v = (y,t,L) G ipi of 
u with R(v ) > 0. We show R(u) + s > R{v) + t by assuming on the contrary that 
R(u) + s < R(v) + t. There exist a k G {0, . . . , i — 1} and a w with u G ^fc+i \ V'fc an d 
{v,w} = Vi{v). The algorithm implies that R(v) + s < d(v,w). Lemma 12.21 and our 
assumption R(u) + s < R(v) + t yield R{u) + s > d(u,v). This gives d(u,v) < d(v,w) 
which contradicts the mutual nearest neighbour property of v and w in ip \ ipk-i- 

In the remaining case u G if) \ ipoo we can proceed as in the second case of the last 
paragraph. Thus every u G if) has at least one earlier neighbour. □ 
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Lemma 2.6. Assume that if) G H. Then R := R(ip, ■) is the unique growth-maximal 
function on if). 

PROOF: Let R' be some growth-maximal hard-core function on if). First we prove by 
by induction over the steps of the algorithm that R(u) = R'(u) whenever u G ipoo satisfies 
R(u) > 0. 

Let u = (x, s, K) G if>x satisfy R(u) > 0. The then yields algorithm that 

d(u, v) = R(u) + s = R(v) + t, (2.16) 

where v = (y, t, L) is the mutually nearest neighbour of u in if). Assume R'(u) = 0. Then 
there exists w = (z, r, M) Eif)\ {u} with R'{w) > a(z, M, x, 0) and s — r > a(z, M, x, 0). 
Hence 

d(u,w) = a(z,M,x,0) +r < s. (2.17) 

If s < t, we get d(u, v) = t + a(x + (t — s)K, K,y, L) > s. If t < s, we also get 
d(u, v ) = s + a(x, K,y + (s — t)L, L) > s. In view of (12.171) we obtain that w ^ v and 
d(u,w) < d(u,v). But w 7^ v leads to d(u,w) < d(u,v). As this contradicts the mutual 
nearest neighbour property of u and v, we must have that R'(u) > 0. 

Assume that R(u) ^ R'{u). By Lemma [2721 and (12.161) we can assume w.l.o.g. that 



R'{u) +s < d(u,v). (2.18) 

Let w G ip be a stopping neighbour of u w.r.t. R'. In particular, R'(w) > 0. By Lemma 
12.21 we conclude R'{u) + s > d(u,w). Hence (12. 18[) implies that d(u,w) < d(u,v), a 
contradiction. 

Now we take u = (x, s, K) G ^j+i \ V'i f° r some i > 1 and assume that R(u) > 0. 
Again we first show that R'{u) > 0. Assume R'(u) = 0. Then there exists w G ip \ {u} 
with R'(w) > a(z, M, x,0) > and s — r > a(z, M, x, 0). Assume w G Then we 

can argue as in the case u G ipi to conclude that R'(u) > 0. Hence we can assume that 
w G ipi. from the induction hypothesis we get R'(w) = R(w). Therefore, 

d(w, u) = a(z, M, x, 0) + r < R(w) + r. (2.19) 

Due to the algorithm < R(w) + r < d(w, w'), where w' is a mutually nearest neighbour 
(w.r.t. R) of w in ip \ ipj for some j < i — 1. Therefore d(w,u) < d(w,w') contradicting 
the choice of w' G if) \ ipj D ip \ tpi and the fact that u G if>\if>i- Thus -R'(m) > 0. 
Next we show that R{u) = R'{u). Assume Vi{u) = {u,v} and 

R(u) = s[(u) — s = d(u, v) — s. 

Assume moreover that R'(u) < s'^u) — s. Then we can argue as at f 1 2 . 1 8 [) to conclude 
that u cannot have an earlier neighbour in if)\if>i w.r.t. R'. Hence there exists an earlier 
neighbour w G if>i of u w.r.t. R', that is R'(w) + r < R'(u) + s and R'(w) > 0. By induction 
hypothesis R(w) = R'(w). Therefore the algorithm implies 

f{u, w) = inf{r > : (x + rK) n {z + R(w)M) ^ 0} + s = R'{u) + s < s'^u). 
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This leads to 



f{uM<s' l {u)<s'l{u) 



inf f(u,v) < f(u,w) < 4(u) 



a contradiction. Now assume R'(u) > — s. The hard-core property of R' yields that 
R'(v) < Si(v) — t. But we have just proved that this is not possible. Hence only the 
desired case R'{u) = s'^u) — s = R(u) remains. 

Now assume R(u) + s = s"(u) < s'^u). Then there exists v G tpi with R{v ) > and 
s"(u) = f(u, v) + s. By induction hypothesis R(v) = R'(v). Hence the hard-core property 
of R' that 



To show that R(u) = R'(u) we assume on the contrary that the above inequality is strict. 
Since, by induction hypothesis, R(w) = R'(w) for all w G ipi with R(w) > 0, we obtain 
from s"(u) + s = f(u, v) and the definition of f(u, v) that u has no neighbour in ipi w.r.t. 
R'. Because of s"(u) < we get that < R'{u) + s < f(u,v) < d(u,v') where v' is 
the mutual nearest neighbour of u in ipi+i \ ipi. By Lemma [2.21 we deduce that 



for every neighbour w = (z, r, M) G ip \ ipi of u. Therefore U cannot have an earlier 
neighbour in ip \ ipi w.r.t. R'. This is a contradiction. 

Next we show R(u) = R'(u) for u G ip \ ipoo- Due to the algorithm we can assume 
R(u) > 0. Because of R(v) = R'{y) for all v G ipoo with R(v) > and the hard-core 
property of R! we obtain the inequality R'(u) < inf t)6 ^, oo) ii( t) ) > Q f(u, v) — s = R{u) . If 
R'(u) < R(u) the point u would not have an earlier neighbour by the definition of f(u, v). 

Finally we show that R'(u) = holds for all u G ip with R{u) = 0. For u G ip with 
R(u) = let v G ip be an earlier neighbour of u in ip w.r.t. R. Therefore it follows 
R{v) > and x G R(v)L. We have already shown above R'iy) = R(v). Thus x G R'(v)L. 
Since R 1 is a hard-core function this yields R'(u) = 0. □ 

For a Borel set AcK d Vd(A) the volume of A. Further let Kd := Vd(B d ) denote the 
volume of the unit ball. The proof of Lemma [2.41 uses the following lemma. 

Lemma 2.7. Letr > 0, < s < t andK x ,K 2 G K d withK x C x x J rr x B A , K 2 C x 2 + r 2 B d 
for some r 1; r 2 > and Xi,x 2 G M d . Then 



s';(u) + s = f(u,v)>R'(u). 



R'(w) + r < d(u, w) > d(u, v') > R'(u) + s, 



ViirK-L + tK 2 ) - V d (rK x + sK 2 ) < max (f 




PROOF: Using mixed volumes (see e.g. [T2j (5.1.26)]), we obtain 




where we have also used the homogeneity part of [12, (5.1.24)]. Therefore 



V d {rK x + tK 2 ) - V d {rK x + sK 2 ) = J2( ■) r j (t d ~ j - s^ViK^j], K 2 [d - j)). 



j=0 
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Since Ki C x-i + rjl? d we can now use the monotonicity property [121 (5.1.23)] and the 
translation invariance of mixed volumes to obtain that 

d-1 / a 

^(riTi + tAT 2 ) - V d [rK x + s# 2 ) < £ f . V(^' - s^Vi^V^S*), 

i=o 

where we have used that V^(.B d , . . . , 5 d ) = Vd(B d ). This implies the assertion. □ 

Proof of Lemma 12.41 We construct a measurable set G G H for which we can show 
that P(\l/ G G) = 1. Consider an increasing sequence hi > 0, i G No, with lim^oo b{ = oo 
to be specified later and define Wj := [— j,j] d , j > 1. Moreover for n G N, s < i 
and a bounded Borel set S C M d the set F n (s,t,B) is defined as follows. A set ip G 
N belongs to F n (s,t,B), if it contains n + 1 different points tt , . . . , w„ such that s < 
d(uo, u\), . . . , d(u n _i, w n ) < t and u G 5' := B x M + x /C d . Define 

F(s,t,5) := f|F n ( Sj t,S). (2.20) 

n>l 

and G as follows. A set ip G N belongs to G if it satisfies (12.71) and 

card{f E.ip\ {u} : d(u,v) < r} < oo, r > 0, u G ip, (2-21) 
a(x,K,y,L)^\s-t\, (x, s, K), (y, t, L) G ip, x ^ y, (2.22) 

as well as ip ^ F where 

oo oo 

F :=U(J F^b^Wj). (2.23) 

i=0j=l 

We recall the definition of the set H from the proof of Proposition 12.31 as the set of all 
ip G N which satisfy (12.71) and the algorithm of Proposition 12.31 yields (based on ip) a set 
tpoo with card(?/> \ ^oo) < 1 and {v G ipoo : R(ip,v) > 0} ^ 0. We now show G C H 
an take ip E G. Define (p :— ip \ ip^ and assume card<£> > 2. Furthermore we take 
u G <£> and define z/(w) as the set of all mutual nearest neighbours of u in and set 
u{u) = {v}. Define A(u) := {w G ?/> : d{u,w) < d(u, v)} U {w E ip : d(v,w) < d(v,u)}. 
By construction A(u) contains no points of ip \ {v{u) U {u}). Because of (12.211) the set 
A(u) is finite. So it contains no points of ipoo, since the corresponding growth times 
would have been determined until the i-th step of the algorithm for some i G No- Since 
a pair of mutual nearest neighbours u and v in p would also be a pair of mutual nearest 
neighbours in ip \ipi at least one point of (p would have attained a growth time, which is a 
contradiction to the definition of <p. In conclusion v(u) has to be empty. By induction and 
the symmetry of d we show that in this case ip must have a descending chain, i.e. an infinite 
sequence U\, u 2 , ■ ■ ■ of pairwise different points of ip with d{u n , u n+ \) < d(u n -i, u n ) for any 
n > 2 (see also [I] and [2]). Since <p contains no pair of mutual nearest neighbours there 
exist three different points Ui,u 2 ,u 3 with d(ui,u 2 ) > <i(u 2 ,-u 3 ) and m 3 being the nearest 
neighbour of u 2 and u 2 the nearest neighbour of u\. Inductively suppose u\, . . . ,u n , 
n > 2, to be n different points in p with u\ being a nearest neighbour of Ui-i satisfying 
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d(ui, u u -i) < d(ui_i,«i_2), i = 3, ...,n. Let u n+1 be a nearest neighbour of u n . This 
implies c?(w n+ i,« n ) < d(u n ,u n ~i) since otherwise, by the symmetry of d, we would get a 
contradiction to the nearest neighbour property. To show it„_i ^ {iti, . . . , w n } we assume 
u n+ i = Uj for some j G {1, . . . , n — 1}. If w n would not be a nearest neighbour of Uj we 
get for some u G ^ 

d(uj,u) < d(uj,u n ) = d(u n+ i,u n ) < d(uj + i,Uj) = d(iij } Uj + i). 

Since ttj+i is a nearest neighbour of Uj this is a contradiction and u n has to be a nearest 
neighbour of But this leads to u n+ i and u n being mutual nearest neighbours contrary 
to our assumption. Hence ip would contain a descending chain and so if> would also contain 
a descending chain. This implies if) G F (defined in (12.23P ) and therefore the contradiction 

Now assume if) G G and {v G ipoo '■ R(v) > 0} = 0. By our construction algorithm this 
is only possible if there exists a sequence of points it, = (xi,ti,Ki) G ipco, i £ No, such 
that ti < d(ui,u i+ i) < t i+1 and £j > for i G N . This leads to if) G F and again to 
^ G G c . 

We show P(* G F c ) = 1 and that the conditions (fUD, (gUTJ) and f[2T22l) have full 
measure, that is P(* G G) = 1. We start by calculating P(\& G F n (s,t,B)) for rt G N, 
s < t and 5 C M. d . Positive constants (depending only on the dimension) occuring in this 
calculation are denoted by Cj, i > 1. Using assumption (I2.2p and recalling the definition 
B' = B x R + x K d we obtain 

F(V eF n (s,t,B)) (2.24) 

< C? +1 y 1{U G £?'}1{ S < d(Mo, Ml), ... , dK-1, «n) < ® Q) n+1 (d(u , • • • , O). 

With /(tto, • • • , w n _i) := 1{mo £ < ^( M o, ^i), . . . , d(u n -2, u n -i) < we conclude 

P(*eFn(M,5)) 

<C? +1 y /K • • • , «„-l)l{s < rfK_ b « n ) < t}(A d ® Q)" +1 (C/( M0 , • • • , «n)) 

<c™ +1 y /(n , . . . , «n-i)l{s < d(u„_i, w n ) < t, t n _i < t„}(A d <g> Q) n+1 (rf( Mo , • • • , Un)) 

+ C™ +1 J f(uo, . . . , U n _i)l{s < M„) < *, t n _x > t„}(A d ® Q) n+1 (rf( Mo , ■ ■ ■ , U n )). 

(2.25) 

We treat the summands separately, omitting the constant c™ +1 and using the definitions 
fjL n : = (A d <g) Q) n , f n := /(it , • • • , m„_i) and (sc*, := u i: i = 0, . . . , n. By (j23J) the 

first integral in (I2.25P is bounded by 



/ n l{s < tn-i + a(ar„_i, #„_i, x n , 0) < t}// n+1 (d(tt , . . . , «„)) (2.26) 

+ y / n l{s < t n + a{x n -i + (t n - t n -i)K n ^K n ^x n ,K n ) < t}fi n+ i(d(u , . . .,u n )). 

(2.27) 
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We recall that Vd(A) is the volume of A C M. d . The integral ( I2.26P equals 

<t,(s- t n+l )+ < a(x n -i, K n _i, x n , 0) < t - t n _i}/x n+ i(d(tt , . . . , u n )) 

= / fnl{t + [(£ - £ n _i)ir n _i \ (S - t n _l) + K n _l]}// n+ l(d(w , • • • , Mn)) 

< (^-^) // /n^(ifn-lK(dK...,«n-l))Q(d(<n,if«)), (2-28) 



where we have used the inequalities (t — c) d — (s — c) d <t d — s d and (£ — s) d <t d — s d for 
< c < s < t. We define 

h(s, t) := max{f - s* : i = 1, . . . , d}, < s < t, (2.29) 

and recall that p(K) is the radius of the smallest all circumscribing K £ /C d . The integral 
f[272T|) equals 

j '/nl{^n-l < ^> ( s — ^n)+ < ^(^n-l + (^n — ^n-l)-^n-l; ^tn ^n) < £ ~~ ^n} 

p n+1 (d(u , . . .,«„)) 
= y / n l{i n _i < t,x n £ x n _i + (i n - t n _i)if„_i + (i - t n )(K„_i - K n ), 

X n £ Xn-1 + (t n ~ t n -l)K n -l +{S~ t n ) + (if n _i - i^ n )}/X n+ i(d(w 0) • • • , ««.))} 
I0n-1 < *}/n [Vrf((t„ - t n ^)K n _ x + (t- t n )(K n ^ - K n )) 

- V d ((t n - t n ^i)K n _ t + (s - t n ) + (X n _i - fQ) p n (d(u , u n - 1 ))Q(d(t n , K n )) 
<jj f n h(s,t)J2*d(^J {t n -t n ^Yp{K n ^)\p{K n ^)+p{K n )) d - 1 



p n (d(u , . . . , w n _ 1 ))Q(d(i n , fQ), (2.30) 

where we have used Lemma 12.71 to obtain the inequality. Now we calculate an upper 
bound for the second summand in (12. 25ft . This can be done as in (12.281) and ( 12.301) with 
the roles of (t n ,K n ) and (t n _i,lf n _i) interchanged in the integrand. So we obtain 



f(uo, • • • , U n _i)l{s < d(Un-i, Un) < t, t n _i > t n }fjL n+1 (d(u , U n )) 

<2{t d -s d ) jj f n V d (K n )p n (d(u ,...,u n ^))Q(d(t n ,K n )) 

+ 2 J J f n h( S ,t)Y,Kd(f) (tn-l-^V^^P^n-O+P^n-l))' 



i=0 

H n (d(u Q , u n ^))Q(d(t n , K n )). (2.31) 
Combining (|2T28|) - (|2T3B with d2T2^1 yields 

P(* £ F n (s, t, B)) < cl +1 h(s, t)c 2 J j / B ((Vi V 1) + (t n V l)) d 

((p(#„-i) V 1) + (p(K n ) V 1)) d fi n (d(u , u n ^))Q(d(t n , K n )). 
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Since the integrand in (I2.24p equals f n+ i, we obtain recursively 

P(tt G F n (s, t, B)) < V d {B)<%cth{ 8 , t) n (2.32) 

„ n 

J l[((U-i V X ) + (*< V V 1) + (P« V l)) d Q" +1 (rf((t , AT ), • • • , (*n, K n ))). 



i=i 



It is convenient to introduce a random vector ((tq, Y ), . . . , {r' n1 Y n )) with distribution 
Q n+1 and to define n := t[ V 1 and p { := p{Yi) V 1 for i = 0, . . . , n. Then inequality (E32D 
can be written as 



n . n 

P(# G F n (s,t,B)) < V d (B)c n 3 c 4 h(s,t) n E[( JJ( Ti _i +Ti)) + 

i=l i=l 

Now we can use the elementary inequality 



dn 



i=l 1=1 i=0 

that holds whenever a , . . . , a n > 1. This gives 
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2d 2d 
Pi 



P(v& g F n (s,t,B)) < V d (B)c^h(s,t) n 2 2nd E[ J]r| 

i=0 

= ^(5)^ C 4^(5,t) n (4 d )-(E[r M p^) 

Assumption (12. 3p implies that 

P(* G F n (s,t,B)) < V d (B)c 5 c n 6 h(s,t) n . 
Now we define 6j := (i/(c& + i G N . By definition fl2T29|) of h 

h(bi- u bi) < 



C 6 + l 



Hence lim n ^ 00 P(^ G F„(&i_i, Wj)) = 0, i, j G N, so that f[2~20l) implies P($ G F) = 0. 

We now show that (12T7D . (l2T2"Tj) and (12T2"2"|) have full measure. The probability that * 
does not satisfy (12. 7p can be bounded by 

E 2j l{d(wo,«i) = ^(mo,m 2 )} + E 2J l{rf(wo, Wi) = d(ii 2 ,'U3)}, 

(2.33) 

where the sum is taken over all triples resp. quadruples of points of $ with pairwise 
different entries. We recall the definition (12. ip of the factorial moment measure v^ k > for 
k G N and that (|2~2|) holds for \1>. Together with the definition (12T5|) of d it follows that 
the first summand of (I2.33j) is bounded by 



a III l{x 2 G dA(u ,ui, (s2,K 2 ))}dx 2 Q(d(s2,K2))Li2{d{uo,Ui)), 
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where A(u , U\, (s 2 , K 2 )) is a convex body, depending on u , U\ and (s 2 , K 2 ) but not on x 2 . 
Hence this integral is zero. The second summand in (I2.33P can be treated analogously. 
Thus we obtain that \I/ satisfies (12 .7p almost surely. That \1/ does also satisfy ( 12 .22 j) almost 
surely can be shown with the same type of argument. 

Finally we prove that P(card{t> G ^ \ {u} : d(u, v ) < r} < oo for all u G = 1 for 
all r > 0. It is sufficient to show, that 

E l{card{v E^\{u}: d(u, v) < r} < oo}l{u G Wj x R + x /C d } (2.34) 

is finite for all j e N and r > 0. For fixed j and r this expectation is bounded by 

l{x G Wj, d((x, s, K), (y, t, L) < r}dxdyQ 2 (d((s, K), (t, L))), 
which in turn is bounded by 

l{x G Wj, s<t,s + a(x, K, y, 0) < r}dxdyQ 2 (d((s, K), (t, L))) 
+a 2 J jj l{x G Wj, s <t,t + a{x + (t- s)K, K, y, L) < r}dxdyQ 2 (d((s, K), (t, L))) 

+fl2 /// e Wj ' t - s,t + L ' x ' °) - r } dxd y&( d (( s ' K )i (*» L ))) 

+a 2 l{x G Wj, t<s,s + a(y + (s- t)L, L, x, K) < r}dxdyQ 2 (d((s, K), (t, L))). 





(2.35) 

The first two summands equal 




l{x G Wj, s < t, s < r}l{y G x + (r - s)K}dxdyQ 2 (d((s, K), (t, L))) 

W Jffn* e W„.<t, t <r}l bl exHt- s)K + (r - t )(K - L)} 

dxdyQ\d((s,K),(t,L))). 

This sum is bounded by 

l{x G Wj, s<t,s< r}K d (r - s) d p(K) d dxdyQ 2 (d((s, K), (t, L))) 




+a 2 JJJ l{x G Wj, s<t,t< r}K d [(t - s)p(K) + (r - t)(p(K) + p(L))] c 

dxdyQ 2 (d((s,K),(t,L))). 

By (12. 3p both summands are finite. By symmetry the third and fourth summand in (12.351) 
can be treated in the same way. Therefore (I2.34p is finite for all j G N and r > 0. □ 
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3 Absence of percolation 



In this section we consider a point process \I/ = {(X n ,T n , Z n ) : n > 1} satisfying the 
moment assumptions (12. 2p . where the probability measure Q is assumed to satisfy (12.31) 
and, moreover, to be concentrated on the set of all strictly convex bodies. In addition we 
assume \l/ is non-empty and stationary, that is F(x + ^ G •) does not depend on x E M d . 
The intensity of ^ is defined by 7* := E*([0, l] d x R + x K d ). By this is a finite 
number, while \1/ 7^ implies 7^ > 0. 

Due to Theorem 12.11 there exists a P-a.s. unique growth- maximal hard-core model 

= {(X n , T n , Z n , R n ):n> 1} based on ^. We define 

00 

Z :=Z{V):=\JX n + R n Z n = (J X n + R n Z n , (3.1) 

n=l n:_R n >0 

that is the union of all grains (which started growing). Note that X n e Z if R n = 0. 
We say that Z (or the growth-maximal hard-core model) percolates, if Z contains an 
unbounded connected component. To be more exact we introduce a graph with vertex 
set ^+ := {(X n ,T n , Z n ) : R n > 0}. Two different points (X m ,T m , Z m ), (X n ,T n , Z n ) G ^+ 
share an edge if they are neighbours, that is if (X m + R m Z m ) fl (X n + R n Z n ) 7^ 0. A 
cluster is a connected component of this graph. 

As P(\& 7^ 0) = 1, the second defining property (ii) of a growth- maximal hard-core 
model given in the introduction together with Theorem 12 . 1 1 implies that almost surely 
is non-empty too. Since \E' + is stationary (by the translation invariance of R) we have in 
fact that P(card\E' + = 00) = 1. Our aim in this section is to verify the following theorem. 

Theorem 3.1. Almost surely there are no infinite clusters. 

This theorem implies in particular that the random set Z does not percolate. Our 
proof is based on some ideas in [I] and [2] . We begin with the following lemma. 

Lemma 3.2. Almost surely any point (X n , T n , Z n ) e has exactly one earlier neighbour. 

PROOF: Assume without loss that \1/ e H with H being the set from Proposition 12.31 
In addition suppose u = (x,s,K) 6 \l/ with R{u) := R(^,u) > where i?(-, •) is the 
function introduced in Proposition 12.31 By growth-maximality u has at least one earlier 
neighbour. If u has more than one earlier neighbour, then for two different neighbours 
v = (y, t, L) and w = (z, r, M) of u one of the following cases must occur: 

R(u) = R(v) + t - s = R(w) + r - s > 0; (3.2) 
R(u) = R(v) +t-s>0, R(u) + s > R(w) + r > 0; (3.3) 
R(u) + s> R(v) + t > 0, R(u) + t > R{w) + r > 0. (3.4) 

We need to show that all three cases (I3.2p -( l3~~4"j) have probability zero. In order to under- 
stand the argument we first illustrate each of the three cases under the additional assump- 
tion t < s < r. Because of the construction algorithm from the proof of Proposition 12.31 it 
follows from (Q that R(u) = a(x + (t- s)K, K, y, L) = a(x, K,z + (r- s)M, M)+r-s. 
In the case (13. 3 j) we have R(u) = a(x + (t — s)K, K,y, L) = a(x,K,z + R(w)M,0). 
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Moreover there exists an n G N such that R(w) can be replaced by an expression de- 
pending on v and a neighbour W\ = (21,7*1, Mi) & {u, v}, where \& n is defined 
in the construction algorithm in the proof of Proposition 12.31 For s > r± we have 
either R(w) = a(x + (s - n)M, M, zi, Mi) or R(w) = a(z,M,zi + R(w 1 )M 1 ,0). If 
s < r\ we have either that i2(w) = a{x,M,z% + (ri — s)Mi,Mi) + (ri — s) or that 
i2(iy) = a(z, M, zi + R(wi)Mi,0). Now R(wi) can be replaced in the same way using 
W\ and a neighbour W2 G \I/ m \ {it, u, wi}, m < n. After a finite number of steps this 
procedure ends and we have a(x, K, z + R(w)M, 0) = a(x, K, z + f(w, Wi, . . . , Wk)M, 0) for 
some > 0. The third case leads to R(u) = a(x, K, y + R(v)L, 0) = a(x, K, z + R(w)M, 0). 
We can apply the same replacing routine as in the second case for R(v) and R(w). 
In the case f)3.2p we have that 

a(x+ (t- s) + K,K,y+ (t- s)+L,L) (3.5) 
= a(x + (r - s)+K, K,z+(s- r)+M, M) + (r - s)+ -(t- s)+ > 0, 

where h + := max{6, 0} for b G R. In the other two cases the preceding replacement 
process yields countable families F n ,n > 1, of functions / : (R d x R + x /C )" — >■ R+ 
such that the following is true. There exist k,l,m > 1, f\ G F^, / 2 G Fi, f 3 G F m and 
(v u ...,v k ) G (^\W) (fc) , (wi...,wO G (^\M) (0 as well as {w[,...,w'J G ^\{u})^ 
such that 

a(x + (t - s)+JC, K,y + (t- s)+L, L) + (t - s)+ = a(x, K,z + fi(vi, v k )M, 0) > 0. 

(3.6) 

or 

a(x, K,y + f 2 ( Wl ,..., wi)L, 0) = a{x, K,z + f 3 {w[, w'JM, 0) > 0, (3.7) 

Using that the set F := U n€ ^F n is countable and assumption (I2.2p on the factorial 
moment measures it is not difficult to see that the probabilities of all three cases (13. 5ft - 
(13. 6p are zero, provided the set of x G R d satisfying fl3.5H - fl3.Tj) is of Lebesgue measure 
zero. In the case A3. Tj) this means that 

J l{a(x, K, K', 0) = a(x, K, K", 0) > 0}dx = 0, (3.8) 

whenever K is strictly convex and the interiors of K' and K" do not intersect. In fact, if 
the previous properties are satisfied and r := a(x, K, K', 0) = a(x, K, K", 0) > for some 
x G R d , then the infimum r = inf{£ > : (x + tK) D (K' U K") ^ 0} is attained in two 
different points in the boundary of K' U K" . It follows from the Lipschitz property of the 
.fT-distance and the strict convexity of K that the set of points x with this property has 
Lebesgue measure 0. The detailed argument (even for sets more general than K' U K") 
can be found in [6]. The cases (13. 5p and A3.6D can be treated similarly. □ 

Two different points (X m ,T m , Z m ), (X n ,T n , Z n ) G ^ + form a doublet if T m + R m = 
T n + R n . This means that the points stop each other mutually. 

Lemma 3.3. Almost surely any cluster contains at most one doublet while any finite 
cluster contains exactly one doublet. 
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Proof: For u = (x,s,K) e f we write R(u) :— R(^,u) and S(u) := s + R(u). 
Consider a cluster 5 C and assume that {u,v}, {u',v'} are two different doublets in 
S. Lemma [3.21 implies that {u, v} fl {u',v'} = 0. (In the following we ignore P-null sets.) 
Assume (without loss) that S(u) > S(u') and let uq, Ui, . . . , u m be a path in such that 
uo = u, u\ {u,v}, u m = u', and v' {u±, . . . ,u m }. Lemma [3721 implies that S(u) < 
S(ui) (otherwise u would have two earlier neighbours v and u\) and then, recursively, 
S(ui) < S(ui+i) for alH G {0, . . . , m — 1}. Therefore S(u m -i) < S(u m ) = S(u') = S(v'), 
contradicting the fact that u' has only one earlier neighbour. Therefore the two doublets 
cannot be connected by a path, so that H can have at most one doublet. If H is finite, 
then any mgH that minimizes the function S on H must be a member of a doublet. □ 

Proof of Theorem 13. It Assume that 5 c $ + is an infinite cluster without a 
doublet. Take u\ G H and let w 2 be the (unique) earlier grain neighbour of «i. By 
assumption we have S'(wi) > 5(n 2 ). Continuing this way we obtain an infinite sequence 
Ui, u%, . . . of different points in \I> such that R{u n ) > and S(u n ) > S(u n+1 ) for all n > 1. 
Now assume that there are s < t such that s < S(u n ) < t for all n > 1. Then Lemma [2T21 
implies that s < d(u n ,u n+ i) < t for all n > 1. We have shown in the proof of Theorem 
12. II that his event has probability zero. 

In view of Lemma 13.31 it remains to prove that there are no infinite clusters with 
exactly one doublet. It is here, where stationarity plays a crucial role. In contrast to [U[2] 
we use the mass-transport principle (see e.g. [8]) 



E / / l{u G B}g( y ^ + ,u,v)ty + (dv)ty + (du) —E l{v G B}g(^ + , u, v)^ +(du)^f +(dv) , 



that holds for all measurable B C M d x R + x K, d and all non-negative measurable functions 
g, provided that g is translation-invariant. For v G define C(i>) = C(\I/ + ,f) C ^ + as 
the cluster containing v. Let be the set of all v G \P+ such that C(f ) is infinite and has 
a doublet. Define g(^ + ,u, v) := 1 if u G \l/oo, i> G C(tt), and if u is the lexicographically 
smallest point of the doublet of C(v). Otherwise, set g(^ + ,u,v) := 0. With this choice 
of g and B := [0, l} d x R + x /C d the right-hand side of is at most E^ + [0, l} d < 7*. 
The left-hand side vanishes if the intensity 7^ of the stationary point process is zero. 
Otherwise has infinitely many points and the left-hand side is infinite. This shows 
that 700 = 0, as asserted. □ 

4 A central limit theorem 

In this section we assume that ^ is an independently marked Poisson process on M. d with 
intensity 1 and mark space [0, 00) x JC d satisfying P-a.s. the restrictions T = and 



for some fixed c > 1. The assumption T = means that all grains are born at the 
same time (taken as without further restriction of generally), while (14. ip implies that 
all growth times are strictly positive and bounded from below and above independently 





(3.9) 



B d CZ C cB d 



(4.1) 
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of the grain shape. Using stabilization arguments ([HI EI]) we shall prove a central limit 
theorem for the growth times R n , neN. 

Let g be a finite kernel from R + x R d x JCq to R d , where /Cq is the set of all if G fC d 
satisfying B d C if C cB d . We assume g to be translation invariant, that is g(t, x, if, A) = 
g(t,x + y, if, A + y) for all y G R d . In addition we let N denote the set of all ■?/> G N, 
such that for all (x, s, if) G ip we have s = and if G /Cq. We use the abbreviation 
(x,K) := (x,0,K) for all (x,0,K) G V an d i> G N and suppress the birth time in all 
other expressions as well. For all x,y G R d and if, L G /Cg define ip( x > K ) := ip U {(x, if)} 
and ip( x ' K )>(y> L ) ;= ^ u {(x, if), (?/, L)} if such a union lies in N and ^ X ' K > ■= tp resp. 
^{x,K),{ y ,L) ._ ^ otherwise. Set 

p{tfj, x, if, A) := g(R(^ x ' K \ x, if), x, if, A)l {cardWl( „,*) )>1} (4.2) 

for all i/j G No, a; G R d , if G /Cq and measurable A G R d , where i? is the function defined 
in Proposition 12.31 With 

fi(^,A):= p(^,x,K,A), (4.3) 

for tj) G No and measurable A C R d we define a measure nijp){-) := (J,(if), •)• Define the 
observation windows W\ := [—1/2, 1/2]^ and W n := n}l d W\ for n > 2. Furthermore set 

Theorem 4.1. Assume there exist a, (3 > snc/i i/iai £/ie translation invariant kernel g 
in (14.21) satisfies the growth bound 

g{t, x, if, R d ) < at' 3 , t > 0, x G R d , if G ACg. (4.4) 

Moreover suppose that f : W\ — >■ R zs a bounded, almost everywhere continuous function. 
Then the limit a^j := lim^oo Var( J w fdpfy/n exists and as n — > oo 

A=( [ f(n- 1 / d x)^ Wn ,dx)-E\ [ f{n- l ' d x)^ Wn ,dx)\)AN^,a,j). (4.5) 

We will derive this theorem from Theorem 2.1 and 2.2 from |10] starting with showing 
that the growth times R(tp( y ' L \y, L) stabilize. Here we borrow heavily from [U] several 
times. 

For ip G N , y G R d and K G K% define 

D(V, y, L) := inf {a(y, L, x, 0) : (x, if) Etp \ {(y, L)}}, (4.6) 

with a as in ( I2.4j) . If {y,L) G t/> then D(i/),y,L) is an upper bound for the growth time 
R(ip,y,L) of the grain (y,L). 

We recall the function d defined in (12.51) and the constant c > 1 from assumption (14.11) . 
Below we will use that 

^||z - y\\ < d(x, K,y,L)<^\\x-y\\, x, y G R d , if, L G (4.7) 
provided that B d G K G cB d and B d G L G cB d . 
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By a (finite) descending chain in ip G N we mean a finite sequence «o, . . . , u n (n > 1) 
of distinct points of ip for which Wj) > d(-Uj, for alH G {1, . . . , n — 1}. For 

ip G N , y G M. d and L G /Cq a triple (x,K,r) G IR d x K d x (0, oo) belongs to the set 
A(ip,y,L) if there exists a descending chain ito, . . - ,u n in ip( y ' L \ such that «o = (y,L), 
u n = (x,K), d(uo,ui) < D(ip( y ' L \y, L) and r = d(u n -i,u n ). We define 

S^, </, L) := B(y, 2cD^ L \y, L)) U (J B(ar, 2cr), 

(x,_ftT,r)G J 4(i/',J/,L) 

if ^ \ {(y, L)} ^ and S(>, y, L) = M d otherwise and S*(V>, y, L) : = S^, y, L) x 

Lemma 4.2. Let ip, p G N 0; y G M d and L e K d . If ip n S*(0, y,L) = ipf] S*(ip, y, L) 
then S(ip,y,L) = S{(p,y,L). 

Proof: We assume \ 7^ since otherwise the result is trivial. Since 

B((y,2cD(ip^ L \y,L)) C S(ip,y,L) we have D{ip^ L \y,L) = D{p^ L \y,L). 

Take (x, K, r) G A(0, y, L). Thus there exists a descending chain it , . . . , u n G ip( y < L ) 
such that m — (VjL), u n = (x,K), d(u ,Ui) < D(ip( y ' L \ y, L) and r = d(u n ^i,u n ). We 
have u m G S*(ip( y ' L \y, L) for 1 < m < n. Therefore w m G y?. So (a;, K, r) G A((p,y,L) 
and hence A(?/>, y, L) C A((p, y, L). 

Consider (x, K, r) G A(<^, L). So there exists a descending chain %...,« n G 
with u = (y,L), u n = (x,K), d(u ,Ui) < D(p( y ' L \y, L) and r = c?(w n _ 1; u n ). We show 
w m ^ ip for 1 < m < n by induction on m. For m = 1 we have d(u ,ui) < D(<pb> L \y,L). 
Using ( 14. 7p with (y,L) and U\ = (xi,Kx) we obtain \\y— xi\\ < 2cD(p)^ y ' L > , y, L). Therefore 
ui G p fl S*(ip, y, L) and hence U\ G ■0. Now assume our assertion holds for u m , 1 < m < 
k-1. Applying (STTJ) with u m = (x m ,K m ) and u m+1 = (x m+1 ,K m+1 ) yields ||x m -x m+ i|| < 
2cd{u m , u m+ i). So Uk G p fl S"*^, y, L) and thus G ■?/> so that A(<^, w, L) C A(ip, y, L). 
Together with the inclusion proved above we deduce A(ip,y,L) = A{p,y,L) and the 
assertion follows. □ 

We recall from the proof of Lemma \2 .41 the definition of a measurable set G C H such 
that P(*gG) = 1. 

Lemma 4.3. Suppose ip G No fl G and (y, L) G K x /Cq. If S(ip,y,L) is bounded, then 
it satisfies 

R(^ L \y, L) = R{ip^ n S*(iP, y, L) U <p, y, L), 
for all finite p G N with p cR d \ S(ip, y, L) x K,' 1 . 

Proof: Suppose S(ip, y, L) is bounded and ip G N . Since S(ip^ y ' L \y, L) = S(ip, y, L) 
we assume w.l.o.g. (y, L) G ip. Moreover because of Lemma 14.21 it suffices to prove the 
equality in the case tp C S*(ip,y, L). Hence we assume this too. We note that since <p is 
finite the set ip' := ip U p is an element of N fl G and define R(x, K) := R(ip, x, K) and 
R'{x,K) := R(ip',x,K). 

Suppose R(y,L) > R'(y,L) and that (xi,Ki) is an earlier neighbour of (y, L) in ip'. 
We assume for now (x\,K\) G ip. Since (xx,Ki) and (y, L) are neighbours we obtain 
(y + R'(y,L)L) fl (x\ + R'(xi, K\)K\) ^ 0. Because of the hard-core property on ip 
and R(y,L) > R'(y,L) it follows R(xi,Ki) < R'(xi,Ki). Suppose (x2,K 2 ) is an earlier 
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neighbour of (x\,K-\) in ip. The hard-core property on ip 1 yields R(x 2 ,K 2 ) > R'(x 2 ,K 2 ). 
Let (x 3 , K 3 ) be an earlier neighbour of (x 2 , K 2 ) in ip' and assume (x 3 , K 3 ) G ip. Continuing 
this procedure leads to a sequence of points satisfying 

R(y, L) > R'(y, L) > R'(x u K x ) > R{x 1: K,) > R(x 2 , K 2 ) > R'(x 2 , K 2 ) > R'(x 3 , K 3 ) . . . , 

(4.8) 

terminating at (x n , K n ) if (x n , K n ) G <p (so n must be odd). The strict inequalities ensure 
that the points of this (possibly) terminating sequence are all different. By Lemma 12.21 
we get 

R(xi, Ki) > d((xi, Ki), (x i+1 , K i+1 )) > R(x i+1 , K l+1 ) (4.9) 

for odd i G N and the same inequalities if we replace R by R' and i G N is even. 
Furthermore R'(y,L) > d((y,L), (x\,Ki)) > R'(xi,Ki) holds. Therefore (14. 8 p implies 

d{{y, L), (x u Kt)) > d(( Xl , K,), {x 2 , K 2 )) > d{{x 2 , K 2 ), {x 3 , K 3 )) >.... 

Thus the sequence forms a descending chain with d((y,L), (xi,K{)) < Dfy^'^, y, L). If 
the sequence terminates at some point (x n ,K n ) G ip' then (x n ,K n ) G S*(ip,y,L) contra- 
dicting ip n S*(ip,y,L) = 0. If the sequence does not terminate then S(ip,y,L) cannot 
be bounded which is also a contradiction. Hence if R(y, L) > R'(y, L) the assertion is 
shown. The case R'(y, L) > R(y, L) is proven similarly, this time starting with an earlier 
neighbour (xi,Ki) of (y,L) in ip. □ 

Let B(y, r) := {x G M d : \\x — y\\ < r} be the ball around y G M d of radius r > and 
B*(y, r) := B(y, r) x K,q. Furthermore define for ip G N H G, y G M d and K G /Cq 

C/(^, y, L) := inf {r > : S(iP, y, L) C B(y, r)} (4.10) 

and U(ip,y,L) := oo otherwise. Because of the Lemma [4.31 we have 

#(V>' L) , y, L) = R{iP^ n £*(y, 2/, L) U y, L), (4.11) 

for all finite <p G No with ip C M d \ B(y,U(ip,y, L) x /Qj. Hence we call {7 a radius of 
stabilization of R at (y, L) with respect to ip. In [ID] the radius of stabilization T of p with 
respect to ^ and K is defined as 

T(ip,K) := inf{r > :p([ip fl -B*(0, t) U(p],0, K, R d ) = p([ip n 5*(0, £)], 0, K, M. d ) 

for all finite ip G N with ^ C (M d \ 5(0, t)) x K%}. (4.12) 

By definition p(ip,0, K,~&L d ) depends on ip only via the radius R(ip(°' K \ 0, K). Because 
of (I4.1ip it follows that T(ip,K) is bounded by U(ip, 0, K). Now we deal with the tail 
behaviour of U($, 0, Z ). 

Lemma 4.4. Let < 7 < 1. Then there exist constants a, (3 > such that 

W>(U(V,0,Z ) > t) < aexp(-/3t 7 ), t > 0. 

T/ie same /io/ds /or [/(v]/^' 21 ), 0, Zo) /or all x G M d , where Z\ is an independent copy of 
Z Q and also independent of $ . 
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Proof: We start with the first assertion. W.l.o.g we assume t > 1. Define U :- 
0, Z ) and D := D(^, 0, Z ). Let < e < l/(d +1). It holds 

t £ l r „ t e 



{ £/> ( }c{/)>-}u {/)<-,[/> 4 



where c > 1 is the constant from the shape condition (14. ip . Define as the set of all 
ip G N , which contain a descending chain (0, K ), U\, . . . , u n such that 

d((0, K ), Ul ) < D(V, 0, tf„) < t7(4c) (4.13) 

and n > t 1 - £ /2. For if) G N and G /Qj such that V {0 '^ o) e N , 0, #„) > i 
and (I4.13P holds we have ip( > Ko ') g -E. To see this, assume on the contrary that every 
descending chain of ip(°' K °^ which satisfies ( I4.13P would consist of n points with n < t 1 ^ 6 /2. 
Since 

0, K ) C B(0, 2c(n + l)D(if>, 0, K )) 

this would imply 

, t l-e 



[/(V, 0, tf ) < 2c(— + 0, tf„) 



by the definition ( OUj) of C/. Since 0, #„) < t e /(4c), t > 1, and < e < l/(d + 1) 
this would lead to the contradiction U(if), 0, K ) < (3/4)i. So we deduce 

> t) < F(D > + p(V°' Zo) g . (4.14) 

Define k := [(ci + l)t £ 7(4c) d ] and Sj := (i/(ci + l)) 1/(i , i > 0, where ci := Vd(S(0, 2c)) and 
["•] is the ceiling function. Hence > t e /(4c). Set m := |~t 1_£ /(2/c)] and n := |~t 1_E /2~|. 
Note that if t — > oo then m — > oo as well. W.l.o.g assume t > 1 so large that m > 3. For 
^ G E take a descending chain (0, Kq), u±, . . . , u n with n > t /2 and thus n > hq. Since 
n > |V~ £ /2] there exists j G {1, . . . , k} and I G {0, . . . , n — m} such that at least m 
consecutive of the n distances d((0, K ), u\) and d(itj_i, Ui), 2 < i < n, lie in the interval 
Ij := [sj-i,Sj]. In addition 

, „ ft l - £ \2ct £ 

Fill < + 1J — < *, 1 < i < w , 

where is the projection of u% onto the first coordinate, that is U{ = {xi,Kj). Hence 
E C Uj =1 Ej where ip belongs to Ej if it contains a descending chain (xi, Ki), . . . , (x m , K m ) 
with ||xi|| < t and lfi_i), (a;*, Kj)) G i G {1, . . . , m — 1}. Hence 

k 

F(¥°' Zo) e E) <^2f{^ e Ej). (4.15) 

3=1 

Furthermore we get 

P($ GEj) <E[ ^ l{||xi|| < t,d((xi,Ki), (xi+i,Ki+i)) G /j} 

((xi,Xi),...,(x m ,i^ m ))g*('™) 

= J ' ■■■ J l{\\xi\\ < t,d((xi,Ki), (x i+ i,K i+1 )) G Ij}dxi . . . dx m Q(dKi) . ..Q(dK m ). 
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A calculation similar to the one following (I2.24p yields 

P(* G Ej) < c cT-H d {s^ - sti)™- 1 = Cb^f— , 

\C\ + 1, 

with a constant cq > 0. Hence there exist constants c 2 , C3 > such that 

P(* e £<,•) < c 3 t d exp(-c 2 m), j = 1, . . . , fc. (4.16) 
Since t ed > 1 and c> 1 we get k < (ci + 2)t £d . Using and ( 147L61) we obtain 

p( 1 jr(o^o) E E)< (c 2 + 2)f d c 3 t d exp(-c 2 m) < c 4 t d+£d exp ( - cgt 1 "^^) 

with constants c 4 , c 5 > 0. Note that 1 — e — ed > 0. Now we give a bound for the first 
summand in (I4.14p . By definition (14.61) of D = D(ty, 0, Z ) we have 

P (d > Q < P(»(fl(o, «^) x «-) = 0) < p(* (*(».£) x Q) = 0), 

(4.17) 

where G /C d , is the largest r > such that 5(0, r) C K . The second inequality 

is due to the first inclusion in (14. ip . Since \1> is a Poisson process this yields 

P [D > \- ] < exp(-c 6 t £d ), 



4c, 

where c@ > 0. With (I4.14p we get from the preceding two inequalities 

F(U >t)< exp(-c 6 t £d ) + c 4 t d+£d exp ( - Cst 1 "^), 

for all sufficiently large t. This implies the first assertion. 

The second assertion can be shown using the same argumentation as above. But to be 
able to deal with the inserted point (x, Z ) we have to divide the event Ej for j = 1, . . . , k 
into events Ejj,i = 1, . . . ,m, where Ejj is the event that Ej occurs and the i-th point 
of the descending chain of the configuration is (x, K) for a fixed K G K d . Hence the 
analogue of (14.151) consists of a double sum on the right-hand side. Since the arguments 
used to bound the probability of Ej also hold for Ej i and because of the definition of in 
we still obtain an sub-exponential decaying tail. We omit the details. □ 

We define 

Pnty, x, K, A) := p(n lA V, n 1/d x, K, n l ^ d A)l Wl (x), 

t(s) := sup F(T(n 1/d ((n- 1/d ^) Wl - x),x,Z Q ) > s), s> 0, 

n>l, x£Wi 

fx p n (A):= Pn((n- 1/d ^) Wl ,x,K,A) 

(x,K)e(n- 1 / d -f) Wl 

for tj; G N , x G M. d , K G JCg and A C R d . Note that 



/ f{n-V d x)n(* Wn ,dx)± [ fdf£. (4.18) 



'Wi 

This equation links (14. 5 p with the form in which the central limit theorem is formulated 
in Theorem 2.2 in [10]. The next lemma shows that (2.7), (2.8) and (2.10) in [10] are 
satisfied. 
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Lemma 4.5. Suppose Z\ is an independent copy of Z and also independent of ^ . Then 



P(E/(tf , 0, Z ) < oo) = F(U(¥ x ' Zl \ 0, Z Q ) < oo) = 1. 
Furthermore for p > 



sup E 

n>l, xeWi 



< OO, 



sup E 

n>l, x,y£Wi 



p„([(n- 1 /^) Wl ] ( ^ l) ,x,Z ,R d ) J 



< oo, 



(4.19) 

(4.20) 
(4.21) 



and sup s>1 s q r(s) < oo for some q > d(150 + 6/p). 

PROOF: The first assertion (14.191) follows directly from Lemma |4~4"1 For p > we have 
sup E[p n ({n- l / d *) Wl ,y,Z ,R d y) = sup E[p(n^-V rf *) Wl X /rf y, £0,^)1 

n>l, j/SWi n>l, j/6Wi 

= sup E[p(^ Wril x,Z ,R d ) p }. 

n>l, X&W n 

The growth bound (14.41) and the definition of p imply for all n > 1 and x G W n that 



E[p(* Wn ,x,Z , W n f] < a p E 



U' 



\(aJ,Z ) 



X,Zo)^l 



{card((* Wn )(^. z o))>l} 



Define the random variable F n |X := R({^ 
we conclude from (2.10) in [9] 



ii'„ 



Zo)l{card(<E< Wn )>o}- For t < diam(W„ 



ny n ,x >t)< p(i?((* 



U' 



\(^o) 



a?,^o) > *) < exp(-c t d ), 



where Cq only depends on W\. Trivially this inequality holds for t > diam(W n ) as well. 
This yields (14.201) . Assertion (14.211) can be shown analogously, if we replace Y n x by 

y w := R((* w J x ' z oUv< z i\x,Z ). 

To prove the statement about r we note 

t(s) = sup F(T(n 1/d (n- 1/d ^) Wl - n 1/d y, Z ) > s) 



n>l, yeWi 

■■ sup F{T{^ Wn 

n>l, xeW n 



X, Zn) > s) 



< sup F(U(* Wn -x,0,Z o )>s), 

n>l, x€W n 

where we use T((p, K) < U(tp, 0, K) for tp G N , .K" G 

With some minor technical modifications we can use the proof of Lemma 14.41 to show 
that F(U(^w„ — x,0,Zq) > s) is sub-exponentially decaying independent of n > 1 and 
x G W„. The dependence on n vanishes because it is possible to refer to \I/ instead of 
^w„ for all n G N if necessary. The dependence on x vanishes because of the translation 
invariance of the Lebesgue measure. We omit the details. □ 

Proof of THEOREM l4.lt Because of Lemma l4.5l the assumptions of Theorem 2.1 and 
2.2 in [10] are satisfied. Due to (I4.18P Theorem 14. II is implied by these two theorems. □ 
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Remark 4.6. The first inclusion in the assumption (14. ip has been used to achieve the 
bound in ( 14.171) . As the proof shows, this bound can be replaced by the slightly weaker 
assumption 

F(tB d C Z Q ) > 1 - exp(-(l/t) £ ) 
for all sufficiently small t > 0. 

Consider a measurable function h : 1R + x JC d — > R+. If we choose g(t,x,K,A) = 
h(t, K)5 X {A) for t > 0, x G M d , if G /Q|, and measurable A C l d we get 

/i(^A) = l{card(^J>2} ^ ^( J R(^,x ) K) ) K)5 x (A)) ) ^ G N . 

Recall that V^(-) denotes the volume and assume, for instance, that h(t, K) = Vd(tK) and 
a constant density / = 1. Then we can use the hard-core property to deduce 

/ f(n-^ d x)ti^f Wn ,dx) = V d (Z(* Wn ))l{caid(* Wn ) > 2}, 

where Z{^>w n ) is the union of the grains of the growth-maximal hard-core model based 
on ^w„ as defined at (13.1 1) . Similar considerations lead to a central limit statement for 
the surface area in the case of strictly convex particles. 
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